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PREFACE 


In computational problems it is common to find that certain functions (such as 
square roots, logarithms, sines, cosines, tangents, etc.) occur either repetitiously in a single 
problem or frequently from problem to problem. These functions are basic functions. 
Programs for computing these basic functions are generally referred to as ‘‘subroutines.”’ 
Because of the frequent occurrence of these basic functions, it becomes expeditious to have 
available complete programs of subroutines and, in some cases, even to prepare paper 
templates of the programs. 


This booklet consists of two sections. 


Section I discusses the uses of subroutines in particular problems and presents the 
techniques employed in programming them and in extending their range and/or accuracy. 


Section II is the main part of the booklet, the chief purpose of which is to con- 
solidate in handy reference form some of the more useful subroutines for the Burroughs 
E101 Electronic Computer. The section contains the subroutines employed in evaluating 
the basic, more useful functions and gives programmed examples of techniques which 
increase the utility of the E101. 


At the end is added, for further reference, a short bibliography of handbooks con- 
taining the series for common functions and of sources for general studies in functional 
approximation. 
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SECTION I 


Uses of, and Techniques for, Subroutines 


A. Uses of Subroutines 


Some of the various situations which arise in the handling 
of subroutines and some programming methods which can 
be considered for these situations are: single use of a sub- 
routine; multiple use of a subroutine; constant increments; 
compound functions; and multi-purpose use of subroutines. 

Before control is transferred to the subroutine, one must 
be sure that the “starting point” requirements given in the 
pinboard program for the subroutine have been met. 


1. Single Use of a Subroutine 


When a subroutine is used only once in the solution 
of a problem, there are two procedures, either of which 
may be followed: 


(a) to incorporate the subroutine as part of the main 
program (the sequence of steps in the subroutine 
would not be altered, although the steps need not be 
located on a single pinboard; the constant required in 
the subroutine may be stored in any convenient mem- 
ory addresses) ; 


(b) to use a separate pinboard for the subroutine and 
one “U”’ instruction to transfer from the main part of 
the problem to the subroutine pinboard and a second 
“U” instruction to transfer back again. 


The first of these two procedures fits the subroutine 
into the main part of the problem and, as a result, econ- 
omizes programming space and running time. The second 
is a bit more flexible and makes better use of pinboards 
or templates on hand. 


2. Multiple Use of a Subroutine 


The multiple use of a subroutine is advisable when 
a subroutine is to be used in several places in the solu- 
tion of a problem. Each use of the subroutine involves 
transfers, first, from the problem to the subroutine and, 
then, back to a different part of the problem. The E101, 
with the aid of a “subroutine exit instruction,’ has a 
simple procedure for handling these transfers. 

Each call for the subroutine requires two steps in 
the main program: 


(a) “H 0b.” (This instruction homes the E switch to 
“b,” the “‘b” being pinned to correspond to the pin- 
board to which the program should return after the 


use of the subroutine) ; 

(b) “Ua 0.” (This instruction transfers control to the 
beginning of the subroutine, the ‘“‘a’’ being pinned to 
correspond to the pinboard which contains the sub- 


routine). 


Step No. 


Neither of these two instructions affects the contents 
of the accumulator. 

The transfer back to the appropriate step in the main 
program, i.e., the step following the “U a 0” instruction, 
is effected by the “subroutine exit instruction,” “U E *.” 
The “U E *” instruction, which occurs as the last step in 
the subroutine, transfers control back to pinboard “‘b” 
of the last ‘‘H 0 b” instruction of the main program and 
then to the step following the ‘U a 0” instruction on the 
same pinboard “b.” 

To illustrate the procedure, let us assume that we 
wish to transfer to the sine subroutine (pinboard 4) 
from the main routine, which is first at step 5 of pin- 
board 1, later at step 4 of pinboard 2, and finally at step 
3 of pinboard 3. After each use of the sine subroutine in 
pinboard 4 we must return to the pinboard (1, 2, or 3) 
from which the transfer was made. 

Specifically, in the illustration below, the multiple 
use of the sine subroutine would necessitate initially a 
transfer from pinboard 1 (step 6) to the beginning of 
the sine subroutine on pinboard 4 (step 0); then from 
the end of the sine subroutine on pinboard 4 (step 15) 
back to the main routine on pinboard 1 (step 7); later 
from pinboard 2 (step 5) to the beginning of the sine 
subroutine on pinboard 4 (step 0); then from the end 
of the sine subroutine on pinboard 4 (step 15) back to 
the main routine on pinboard 2 (step 6); later from pin- 
board 3 (step 4) to the beginning of the sine subroutine 
on pinboard 4 (step 0); then from the end of the sine 
subroutine on pinboard 4 (step 15) back to the main 
routine on pinboard 3 (step 5). 


Pinboard No. 


subroutine 
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3. Constant Increments 


When it is desirable to find a function of a variable 
that increases at a constant rate, the “chain formula” 
method should be considered. As an example, consider 
sin ku, and cosku, (used frequently in Fourier Series 
and in other applications, k being an integer that steps 
from 1 to n). One method would be to use, by means 
of two pinboards, the sine and cosine subroutines, shown 
on pages 6 to 8. 

A better method, however, is the constant-increments 
one, which makes use of the chain formulas for sin ku, 
and cos kuy: 


sin kuy = sin [(k — 1) + 1]u, 

= COS Uy sin (k — 1)uy + sin uy cos (k — 1)uy 
and 
coskuy = cos [(k — 1) + 1]up 

= COS Uy cos (k — 1)uy — sin uy sin (k — 1)uy. 


Using these formulas, we can program a subroutine 
that computes both sin ku, and cos ku, in less than one 
pinboard. The program requires 14 steps and 2 seconds 
for each value of k as compared to 26 steps and 8 seconds 
required in the two-pinboard method (see page 18). 

If only cos ku, is desired, the chain formula, 


cos kup = 2 cos (k — 1)upy cos up — cos (k — 2)uUp, 


may be used. 

Another example in which the constant-increments 
approach may be used is es. Once e%» has been com- 
puted, e'. can be found very simply for any number of 
integral increments merely by adding a few steps to the 
program, i.e., by raising e%. to successive integer powers. 

The constant-increments approach can be used effec- 
tively with other functions whenever the argument of 
the function is to be advanced by multiples of an integer. 
Because of the word length of the E101 (12 digits), 
there is rarely, if ever, a build-up of errors sufficient to 
cause trouble in engineering computations. 


4. Compound Functions 


When a problem involves a compound function, 
such as 
e"sinu or sinhu + logu, 


in which each of the component functions can be ex- 
pressed as a power series, there are two possible methods 
of attack. 

One method is to compute each component function 
separately and then find their combined value. 

Another, and in some cases a better, method is, first, 
to combine the power series of each component function 
into a single new power series and, then, to program the 
new series. As an example of this method, we can multi- 


ply the e" series by the sinu series to arrive at a new 
series (see page 19). The program for the e"sinu 
series requires only 8 steps and 10 seconds as compared 
to 26 steps and 16 seconds required in the two-pinboard 
operation. 


5. Multi-Purpose Use of Subroutines 


When a problem involves several functions, each of 
which can be approximated with sufficient accuracy by 
a finite polynomial, an economy in the number of pin- 
board steps required can be effected by the use of a 
multi-purpose subroutine, i.e., one that will suffice for 
all the functions involved. Such a multi-purpose sub- 
routine for sinu, cosu, arctanu, and e~", illustrated on 
pages 19-20, requires only 10 steps as contrasted with 
the need for 4 pinboards to compute the functions sep- 
arately. 

Since all four functions can be expressed as poly- 
nomials in u, with only the coefficients varying, the same 
pinboard steps can be used for the polynomials in u, 
with different memory addresses for the coefficients. 
The coefficients for each function are stored in a separate 
row in the memory, the E switch being used to select the 
proper row. 


. Techniques for Subroutines 


1. Extension of Range and Accuracy 


The programs for the subroutines given in Section II 
are valid, within the accuracy specified, only for the 
indicated range of the argument. Since, however, there 
are some problems that involve ranges of the argument 
greater than those indicated in the subroutines in Section 
II, additional programming can sometimes be used to 
cover the extended range. 

Problems requiring an extension of range and/or an 
increase in accuracy are too varied and numerous, even 
to catalogue. Therefore, only general lines of approach 
can be indicated. 

An example of a subroutine that can be easily ex- 
tended in range is the exponential function, e*. The 
basic program on page 12 is valid for —1 <u < 1, but 
the range can be extended by the use of either of the 
following identities: 

ev = — or et (675). 
Furthermore, greater accuracy of e" can be obtained if, 
for example, e%/* is computed and then raised to the 
fourth power. 

An identical approach can be employed to extend 
the range and/or accuracy of the 10" subroutine. 

Another example of a subroutine that may be easily 
extended in range is the sinu subroutine, where the 
range, —1/2 <u <7/2, may be increased by the use 


of trigonometric identities, which replace u by an angle 
y in the first quadrant: 


sinu = sin (Z sty y) = cosy, 


sinu = sin (7 + y) = — siny, 


sinu = sin es 3F y) = — cosy, 
sin u = sin (27 + y) = siny. 


If u is a large angle, a standard procedure is to write 


sinu = sin 2 ce eee al 
ae Qa |u| Qa 3 


an | 
where the symbolism E 


jul ,, 


means ‘‘the largest integer in Sai 
Tv 


Trigonometric identities can usually be employed to 
extend the range of any subroutine for a trigonometric 
function. For example, the range for the tan u subroutine 
can be significantly extended by use of the following 
identity: 

ais =o : 
1 Stanza 

In any of the subroutines that use Taylor Series, the 
range and/or accuracy may be extended by increasing 
the number of terms in the polynomial approximation. 
For example, to extend the range of the sin u subroutine, 
on page 6, to —7 <u < 7 and retain the accuracy of 
+0.000 05, the following polynomial would be used: 

¢ u? u> u? u? yi uis 
pea Tike Ti ev mo ae 

Many of the subroutines in this booklet use poly- 
nomial approximations not derived from Taylor Series. 
It is nonetheless true that the range and accuracy of these 
subroutines may be extended by increasing the number 
of terms in the polynomial approximation. Hence, in 
order to aid the user in extending the range and/or 
accuracy of these subroutines, references which contain 
extensions of the polynomial approximation are given 
with each subroutine. 


2. Programming 


The subroutines in this booklet are based primarily 
on polynomial evaluation. The polynomials are written in 
terms of u, u2, or (u — a)/(u + a). The general method 
for programming such an approximating polynomial, 


P(u) = S ay Uy, 
k=0 


can best be shown by an example. 


Suppose n = 4. Then the polynomial becomes 
P(u). = 454+ aye asu? + a,u® + amt, 
which may be written as 

P(u) = a + ula, + u[a, + ula, + au]]]. 

This factored expression is the form leading to the most 
efficient programs for digital computers because [a, + 
a,u] involves merely a, times u, added to a;, which result 
becomes the new coefficient of u in [a, + u[a, + a,u] ]. 
The process can be repeated, and we are thereby able to 
evaluate a high-degree polynomial by repeatedly perform- 
ing a simple multiplication, followed by a simple addi- 
tion. 

Now suppose that the a’s and the value of u, for 
which the polynomial is to be evaluated, are stored in 
memory addresses as follows: a, in 20, a, in 21, a, in 22, 
a, in 23, a) in 24, and u in 25; then the original poly- 
nominal P(u) can be easily evaluated by employing its 
factored form in the following program: 


Step Instruction Contents of Accumulator 


0 R20 a, 

i B 

2 PS) a | a,u 

3 +21 a, + au 

4 B 

5 ete date u[a, + a,u] 

6 +22 a, + ula, + a,u] 

wh B 

8 Ku .5 u[a, + u[a, + a,u]] 

9 + 23 a, + ula, + u[a, + a,u]] 
10 B 

ee oe ula, + ula, + ula, + au]]] 
12 +24 a) + ufa, + ula, + ula, + a,u]]] = P(u) 
13 r Print P(u). 


The foregoing program is linear (not looped). To 
save instruction space we can make use of the “H,” “'S,” 
and ‘“U”’ instructions to form a loop and thereby arrive 
at the following shorter program: 


Step Instruction Contents of Accumulator 


0 R20 

f Head at step 4, when 

2 B F=1, isa, 4-a,u; 

3 X25 F= 2, is a, + ula, + a,u]3 

4 +2F F=3, is a, + ula, + ula, +a,u]]; 

5 $14 F=4,isa,-+ufa, + ula, + ula, + a,u]]]. 
6 U.O2 

7 P Print P(u). 


So far, scaling has been neglected, but now we must 
consider the scaling that would be involved in such a 
problem. First, we must establish some?rules of notation. 


Suppose, for example, that we have a number, say 
20, which we want to write in a memory address. If we 
enter this number into the machine as .200 000 000 00, 
where the decimal is the machine decimal, then we say 
that the given number, 20, times 10-?, is in the memory. 
If we enter this number, 20, into the machine as 
.020 000 000 00, where the decimal is the machine 
decimal, then we say that the given number, 20, times 
10-°, is in the memory. In general, if we enter a given 
number into the E101 in such a way that the machine 
decimal is p places to the left of the real decimal, then 
we say that the given number times 10-? is in the mem- 
ory, B-register, or accumulator, as the case may be. 

If, in the evaluation of the polynomial, u times 10 
is in the memory, then a,u and a, must be made, by cor- 
rectly scaling a, and a;, to agree in their power of 10, 
since they are to be added. For example, if a, times 108 
is in the memory, then a, times 10%? must be in the 
memory since a,u times 108? will be in the accumulator 
and must be added to a;. By following this procedure 
we can determine the appropriate power of 10 at which 
each of the a’s should be written. (It is obvious that the 
printed answer will be of the same power of 10 as the 
constant ay.) 


Thus, in our example, if u times 10-1 is in the mem- 
ory and if the a’s have been stored as follows: a, times 
10°, a, times 10-1, a, times 10-2, a, times 10-3, and a, 
times 10-4, then all the steps used in evaluating the poly- 
nomial can be correctly performed. 


Scaled Coefficients Partial Sums 

a, times 10° a,u times 10-1 

a, times 10-1 [a,u + a,] times 10-1 

a, times 10-2 [a,u -++ a,]u times 10-2 

a, times 10-8 [[ayu + a,]u- a,] times 10-2 
a, times 10-4 [[ayu + a,]u + a,]u times 10-3 


[[[ayu + a,]u-+ a,Ju- a,] times 10-3 
[[[a,u + a,Ju-+ a,Ju-+ a,]u times 10+ 
P(u) times 10-4 = 
[[Lagu-+ a,Ju-+a,Ju +-a,Ju-+ a, times 10-4 
It is perhaps most efficient in the programming of 
polynomial evaluation to begin the scaling with the co- 
efficient of the term involving the highest power of u, 
say a,, and then to arrange the scale factors of the other 
coefficients, namely, ay, 4-2, an-3)°**,@, to be the 
same, respectively, as the scale factors of a,u, a,u?, 
ad® cs * ase. 


SECTION II 


The Commonly Used Subroutines 


A. Preliminary Information 


This section contains programs of the more useful sub- 
routines and furnishes programmed examples of 3 tech- 
niques that increase the utility of the E101. 

Approximations by power series or rational functions 
have been used. The maximum error indicated for each 
program is the absolute difference between the function 
and the approximation. 

For each subroutine these items have been given: the 
range of u (i.e., the independent variable of the function 
for which the program is usable), the method by which 
the given function is approximated, the constants used and 
their number, the number of temporary memory addresses, 
the scale factors, the switches used, the time required, and 
the detailed pinboard program. 

The constants required for each program are shown 
exactly as they are stored in the memory, but, if in any 
given problem it should be convenient to use memory 
addresses other than those suggested, the memory addresses 
can be changed, provided the program steps referring to 
them are also changed. 

The exit instruction back to the main program has been 
omitted except in the case of the conditional transfer. 

The constants in the subroutines are scaled according to 
the scaling techniques discussed in Section I. Changes in 
the scaling of these constants should be made with caution. 

Emphasis has been placed on economy of program steps 
rather than on economy of running time. Program steps 
have been saved by the use of “S” and ‘‘U” instructions to 
iterate the basic steps of each subroutine. Another approach, 
that of stretching the program out in a linear fashion, in- 
stead of iterating with “S” and “U,” requires less running 
time but more program steps. The linear program for cos u, 
for instance, requires 18 steps and 3 seconds as against the 
iterative (‘‘S’” and ‘“U’’) requirement of only 12 steps (a 
saving of 6 steps) and 4 seconds (a loss of only 1 second 
in time). 


B. Subroutines 


1, 
2s 


sinu (—7/2 <u < 7/2, win radians) 

sinu (—90° <u < 90° or —7/2 Cu<7/2, uin 
degrees or radians, respectively) 

sin A (for large angles, A, in radians or degrees) 

cosu (—27/2 <u < 7/2, u in radians) 

cosu (—90° < u < 90° or —a/2 <u< 7/2, uin 
degrees or radians, respectively) 

cos A (for large angles, A, in radians or degrees) 

tanu (—7/4<u<7/4 or —1 <u<l, u in 
radians 

arccos u and arcsin u (0 <u < 1, result in radians) 
arccos u and arcsinu (—1/2 <u < 1/2) or 

—YV 2/2 <u< Y 2/2, result in radians) 


arctanu (—1 <u <1, result in radians) 


. arctanu (—1 <u < 1, result in radians) 
. arctanu (0 <u < 999, result in radians) 


eG hee) 
eo ¢0 <cu-< 10) 
102 (0: a-< 1) 


» dogg Go <i0'< 10) 


log,u (1 <u < 10) 
log,u (1 <u < 10) 


. sinhu (—4.5 <u < 45) 
. coshu (—4.5 <u < 45) | 
2 tanh (= 2:< a<s 2) 
Wa 

ery 


Multiplication of complex numbers: 
(a + bi)(c + di) for —.9999 < a,b, c,d < .9999 


. Division of complex numbers: 


(a + bi)/(c + di) for —.9999 < a,b, c,d < .9999 


1. sinu 2. sinu 
(u in radians) (u in radians or degrees) 


Range: —7/2 <u< 7/2 Range: —90° <u < 90° or —27/2 Cu< 7/2 


Accuracy: Maximum error: 0.000 05 Accuracy: Maximum error: +0.000 000 1 


' Method: Reference—ElectroData. See also Hastings, pp. 
Method: Reference—Taylor Series 


138-140. 
; w uw u™ vw? Let x = u/90 or u/7/2 
| SE oe gh cal cae then 


Number of Constants: 5 : . 7 : 
sinu = sin > x = by Guage 
k=0 
Number of Constants: 7 
Switches Used: F (4 to 8) Number of Temporary Memory Addresses: 2 
Switches Used: F (3 to 7) 
Time on E101: 4.0 seconds 


Number of Temporary Memory Addresses: 2 


Time on E101: 4 seconds 


Special Feature: easier to extend in range and accuracy 


than other sine routine on this page. Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory - Memory Address 
Constants: 
Co | 0.000 015 148 42 97 
As Constant Assumed 
Constant Scale Factor «Appears in Memory © Memory Address ms = Oe. a. a . 
Cs —1 0.007 968 967 93 94 
—1/1! 3 0.198 412 698 41— 94 
C3 —1 0.064 596 371 11— 95 
1/5! 1 0.083 333 333 33 95 
C1 —1 0.157 079 631 85 96 
—1/3! —1 0.016 666 666 67— % 
1 —3 0.001 000 000 00 97 
ia ; ie Se in il 1/90 2 1.411 111 111 10 98 (for u In degrees) 
bs : ./ar 1 6.366 197 720 00 98 (for u In radians) 
Pinboard Program: Pinboard Program: 


Starting point: u (in degrees) times 10-? or u (in 
radians) times 10-1 is in the accumulator. u should be 
shifted to meet this requirement. 


Starting point: u times 10-1 is in the accumulator. 
u should be shifted to meet this requirement. 


o wWo9i OG: Bb 
1 B Tex 98 
2 xX 91 2 W9I1 
3-" 2W 02 5 aD 
4 R 98 4 xo 
Se 14 5° We92 
G— OB 6° R97 
| 7 X/92 7 H13 
8 + 9F 8 B 
| eG Dx 
10 U 06 10 + 9F 
ite 3B phil S 16 
12 > aed 12 U 08 
13 Iz 3B 
14 te <9) 
15 15 
Result: Result: 
sin u times 10-¢ is in the accumulator. The result can sin u times 10-1 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. then be shifted to meet the requirements of any problem. 


3. sinA 
(for large angles) 
Range: —10,000° < A < 10,000° or 
—100 radians < A < 100 radians 
Accuracy: Maximum error: +0.000 000 2 
Method: Reference—ElectroData 


8 


sin 27 u = > Con igs 


what Gott ee 
ah ~~ Qa |A| 2a 


Number of Constants: 10 


Number of Temporary Memory Addresses: 2 


Switches Used: F (1 to 10) 
Time on E101: 8 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Ciz a 0.000 606 752 76 90 
Cis =—9 0.006 406 975 47— 91 
Cis = 0.037 414 651 35 92 
C11 =) 0.150 465 962 59— 93 
Co =) 0.420 407 965 50 94 
Cr a! 0.767 019 335 43— 95 
Cs ~2 0.816 047 826 56 9% 
C3 ay 0.413 416 770 70— 97 
C1 = 0.062 831 849 10 98 
1/360 3 2.117 777 777 80  99(for.A in degrees) 
1/2ar 1 1.591 549 430 00 99(for A in radians) 


Pinboard Program: 


Starting point: A (in degrees) times 10-* or A (in 
radians) times 10- is in the B register. A should be 


shifted to meet this requirement. 
99 
12 
21 


88 
88 
89 
90 
11 


89 
oF 
19 
09 
88 


DOO Oe SAA he eS 
XAGetxX etn ax eo > > 


a’ 
WM od Ww 


Result: 


sin A times 10-? is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


4. cosu 


(u in radians) 
Range: —7/2 <u < 7/2 
Accuracy: Maximum error: +0.000 05 


Method: Reference—Taylor Series 


u2 ut u& u& 
CS SE enn ee oa ae 


Number of Constants: 5 

Number of Temporary Memory Addresses: 1 
Switches Used: F (4 to 8) 

Time on E101: 4 seconds 


Special Feature: easier to extend in range and accuracy 
than cos routine on p. 8 


Constants: 
As Constant Assumed 

Constant Scale Factor Appears in Memory Memory Address 
—1/6! 2 0.138 888 888 89— 84 

1/4! 0 0.041 666 666 67 85 
=—172 =) 0.005 000 000 00— 86 

] —4 0.000 100 000 00 87 

1/8! 4 0.248 015 873 02 88 


Pinboard Program: 


Starting point: u times 107 is in the accumulator. u 
should be shifted to meet this requirement. 


Oo W 89 
i B 
288 
3 W 89 
“4 R 88 
5 H 14 
65/25 
Lao 
So aecor 
2 S 5 Av 
10 U 06 
11 

12 

13 

14 

15 


Result: 


cos u times 10-4 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


5. cosu 
(u in radians or degrees) 


Range: —90° <u < 90° or —7/2 Cu<7/2 
Accuracy: Maximum error: +0.000 000 1 


Method: Reference—ElectroData 
Let x = u/90 or u/7/2 


5 
Tw 
Then cos u = cos Si = ) Gy x2* 
k=0 


Number of Constants: 7 

Number of Temporary Memory Addresses: 1 
Switches Used: F (4 to 9) 

Time on E101: 4.5 seconds 


Constants: 


As Constant Assumed 
Constant Scale Factor Appears In Memory Memory Address 


1/90 2 1.111 111 111 10 83 for u in degrees) 


2/ar 1 6.366 197 720 00  83(foru in radians) 
Cs =| 0.000 091 785 85 84 
Ce = 0.002 086 279 50— 85 
C4 | 0.025 366 935 70 86 
Co sa 0.123 370 053 81— 87 
Co =] 0.100 000 000 00 88 
Cio =| 0.000 002 388 30— 89 


Pinboard Program: 
Starting point: u (in degrees) times 10-* or u (in 
radians) times 10-1 is in the accumulator. u should be 
shifted to meet this requirement. 


OF 5B 

Ls 25535 
2 WeS2 
3. 

4 X 82 
5” We82 
6 R 89 
7 %#H 14 
8 B 

9 xX 82 
10: ~=—-SF 
11 S 18 
12, Uc08 
13 

14 

15 


Result: 
cos u times 10-1 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


6. cosA 
(for large angles) 
Range: —10,000° < A < 10,000° or 
—100 radians < A < 100 radians 
Accuracy: Maximum error: +0.000 000 7 
Method: Reference—ElectroData 
8 
cos 27 u = a cys 
k=0 
where u = =, ariel [LAL] 
20 | A | 2a 
Number of Constants: 10 
Number of Temporary Memory Addresses: 1 
Switches Used: F (1 to 9) 
Time on E101: 7 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Cie =#] 0.001 590 159 38 90 
C14 =) 0.015 108 285 .35— 1 
Ci2 ree 0.077 136 624 40 92 
Cio =2 0.263 212 932 60— 93 
Cs 74 0.602 103 183 15 94 
Ce =) 0.854 504 972 05— 95 
C4 —P 0.649 388 111 90 96 
Co =e O97 390-8810 51 97 
Co =2 0.009 999 998 80 98 
1/360 3 2.111 7117 7717 80  99( for A in degrees) 
1/2ar 1 1.591 549 430 00 99(for A in radians) 


Pinboard Program: 
Starting point: A (in degrees) times 10-* or A (in 
radians) times 10-* is in the accumulator. A should be 
shifted to meet this requirement. 


29 
12 
21 


89 
89 
89 
11 
90 


Se 
NF COMO ON AYU RWwWN YS CS 
n+Xwormsex dgur>xw 


Pee 
WU od Ww 
e 
jo) 
_ 
Oo 


Result: 
cos A times 10-2 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


7. tanu 
(u in radians) 


Range: “ —7/4<u< 7/4 
m —1cucl 


Accuracy: Maximum error: “ +0.000 000 7 
(>) +0.000 2 


Method: Reference—Taylor Series 
tan'u =u + cu® + cu’ cu = += + cu? 


Number of Constants: 10 
Number of Temporary Memory Addresses: 2 
Switches Used: F (0 to 8) 


Time on E101: 9 seconds 


Constants: 
As Constant Assumed 

Constant Scale Factor Appears in Memory Memory Address 
Cs 0 0.000 590 027 44 90 
Cz 0 0.001 455 834 38 1 
Ce 0 0.003 592 128 00 92 
Cs 0 0.008 863 235 50 93 
C4 0 0.021 869 488 50 94 
C3 0 0.053 968 253 00 95 
Co 0 0513374333333) tsa 96 
Ci 0 0.333 333 333 33 97 
1 0 1.000 000 000 00 98 
Co 0 0.000 239 129 11 99 


Pinboard Program: 


Starting point: u times 10° is in the accumulator. u 
should be shifted to meet this requirement. 


0 W 89 
a B 

2 P89 
ele ele, 

4 xX 99 
5 “Halo 
6.) heer 
7 W 88 
8 xX 88 
2 S 17 
10) © U:306 
sg ts, 
12 AYE 
13° 355B 

14 X 89 
15 


Result: 
tan u times 10-1 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


8. arccos u and arcsin u 
(result in radians) 
Range: 0<u<1l 
Accuracy: Maximum error in either function: 
+0.000 000 2 
Method: Reference—Hastings, page 162 


6 

arccos u = (> cs) Vi-u 
k=0 

arcsin u = 1/2 — arccosu 

Number of Constants: 9 or 10 

Number of Temporary Memory Addresses: 3 

Switches Used: F (3 to 9) 

Time on E101: 17 seconds for arccos u 


Constants: 
As Constant Assumed 

Constant Scale Factor Appears in Memory Memory Address 

*ar/21 —1 0.157 079 632 68 50 
(necessary only 
for arcsin u) 

1/2 0 0.500 000 000 00 51 

=) —2 0.010 000 000 00— ~ 52 

C5 0 0.011 146 229 40— 53 

C4 0 0.026 899 948 20 54 

C3 0 0.048 802 504 30— 55 

Co 0 0.088 755 628 60 56 

C1 0 0.214 585 264 70— 57 

Co 0 1.570 796 172 80 58 

Ce 0 0.002 295 964 80 59 


Pinboard Program: 


Starting point: u times 10° is in the accumulator. u 
should be shifted to meet this requirement. 


Pinboard 1 Pinboard 2 
0 W 49 x 51 
Tee PRy 59 W 49 
2 Ee 3 A 12 
5 oD B 
4 %X 49 R 49 
5st SE + 47 
6> “S18 x 51 
7 U 03 + 47 
8 W 48 S 18 
9 R 49 U 03 
10 A 22 B 
Me S52 xX 48 (arccos u times 10-*) 
12 A 5 *A 5 
13 H 10 ¥-. 50 
14 +B (arcsin u times 10-*) 
15 U 20 
Result: 


The result (0 < arccosu < 7/2 or O < arcsinu < 
2/2) times 10-1 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


*These items necessary only for arcsin u. 


9. arccos u and arcsin u 
(result in radians) 
Range: “® —1/2 <u<1/2 
) —V¥2/2<u< V2/2 

Accuracy: Maximum error in either function: 

(a) +0.000 000 1 

(») +0.000 3 
Method: Reference—Langdon, page 13 


4 
arcsin u = as Coy gue 
k=0 
arccos u = 7/2 — arcsinu 
Number of Constants: 5 or 6 
Number of Temporary Memory Addresses: 2 
Switches Used: F (0 to 4) 
Time on E101: 5 seconds for arccos u 
Special Feature: much shorter over a limited range than 
routine on p. 9. 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Cz —1 0.003 822 898 56 70 
Cs —1 0.007 574 473 31 i 
C3 —1 0.016 663 218 82 72 
C1 —1 0.100 000 044 17 73 
Co —1 0.005 315 614 62 74 
*ar/2 | 0.157 079 632 68 15 


(necessary only 
for arccos u) 


Pinboard Program: 
Starting point: u times 10° is in the accumulator. u 
should be shifted to meet this requirement. 


0 W 76 

1 2B 

2 ao 

ee NAST 

4 R74 

5 H 10 

6B 

Ts an 

Sime ee 

DS: 515 

10 U 06 

11 B 

12 X 76 (arcsin u times 10-*) 
“13° AD 
*14 + 75 (arccos u times 10-7) 
15 


Result: 

The result (—7/6 < arcsinu < 76 or 7/3 < arc- 
cosu < 27/3) times 10+ is in the accumulator. The 
result can then be shifted to meet the requirements of 
any problem. 


*These items necessary only for arccos u. 


10. arctanu 
(result in radians) 


Range: —1 <u<l 
Accuracy: Maximum error: +0.000 002 


Methods: Reference—Hastings, page 135 


5 


arctan w= ) Cae gees 
k=0 


Number of Constants: 6 
Number of Temporary Memory Addresses: 2 
Switches Used: F (0 to 4) 


Time on E101: 6 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Cr —1 0.011 643 287 00— 90 
Cs —1 0.019 354 346 00 +a 
C3 —1 0.033 262 347 00— 92 
Ci —| 0.099 997 726 00 93 
C11 | 0.001 172 120 00— 94 
Co —1 0.005 265 332 00 95 


Pinboard Program: 


Starting point: u times 10° is in the accumulator. u 
should be shifted to meet this requirement. 


ie 

1 W 96 
2 oe 8 
3, 8 

4 A 
5 3 
6 H 10 
7 or 
8. -X-97 
2 = ae 
10 «-$: 13 
11 U 87 
i2 6s 

13 X 96 
14 

15 


Result: 


arctan u times 10-1 is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 


11. arctan u 
(result in radians) 


Range: —1<u<l 
Accuracy: Maximum error: +0.000 08 


Method: Reference—Hastings, page 133 
3 


arctan u = Cote 


k=0 


Number of Constants: 5 
Number of Temporary Memory Addresses: 3 


Switches Used: F (0 to 4) 


12. arctan vu 
(result in radians) 
Range: 0 <u < 999.0 
Accuracy: Maximum error: +0.000 09 


Method: Reference—Hastings, page 133 


= u — 1\2** 
arctan u = 7/4 + ) can( 254] 
k=0 7 


Number of Constants: 7 

Number of Temporary Memory Addresses: 3 
Switches Used: F (0 to 4) 

Time on E101: 6 seconds 


Special Feature: true for greater range of argument than 
other routine on this page and the one 


on p. 10. 
Time on E101: 5 seconds Peart: 
As Constant Assumed 
Contents: Constant Scale Factor Appears in Memory Memory Address 
is Conta eaaat Cs ail 0.014 627 660 00 80 
Constant Scale Factor Appears in Memory Memory Address G 44 0.032 118 190 00— 81 
Cs —1 0.014 627 660 00 80 C1 —1 0.099 921 500 00 82 
C3 —1 0.032 118 190 00— 81 a/4 —1 0.078 539 800 00 83 
C1 —1 0.099 921 500 00 82 2 —3 0.002 000 000 00 84 
zero 0 0.000 000 000 00 83 1 —3 0.001 000 000 00 85 
C7 —1 0.003 899 290 00— 90 C7 —1 0.003 899 290 00— 90 


Pinboard Program: 


Starting point: u times 10° is in the accumulator. u 
should be shifted to meet this requirement. 


0” “W935 
1 B 

2- = 95 
3 -W 92 
4 W991 
5: , ato 
6 6B 

Te OE, 
8 =. 7Se 
9 “Sh 36 
10 U 06 
11 

12 

13 

14 

15 


Result: 


arctan u times 10-1 is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 


Pinboard Program: 


Starting point: u times 10-° is in the accumulator. u 
should be shifted to meet this requirement. 


OF esp 
ak B 

2 — 84 
Neat gr ee 
4 R93 
Soe ew, 
6X. 93 
7 Wot 
8 W 92 
5S eae (ges) 
10" B 

5 Re oe 
12 + 8F 
13 S 13 
14 U010 
15 


Result: 


arctan u times 10-1 is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 


13. e 


Range: —1 <u<1 
Accuracy: Maximum error: +0.000 005 


Method: Reference—ElectroData 


6 


= Ga* 


k=0 


Number of Constants: 7 
Number of Temporary Memory Addresses: 2 


Switches Used: F (0 to 6) 


Range: 0 <u < 10 
Accuracy: Maximum error: +0.000 03 


Method: Reference—Hastings, page 182 


et = 7p> cot) 


Number of Constants: 10 
Number of Temporary Memory Addresses: 1 
Switches Used: F (0 to 8) 


Time on E101: 8 seconds 


Constants: 
Time on E101: 6 seconds As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Ea Cr 44 0.012 158 000 00 70 
As Constant Assumed Ce +3 0.022 049 330 00 11 
Constant Scale Factor Appears in Memory Memory Address Cs +) 0.027 728 680 00 12 
Cs 1 0.086 805 555 50 60 Ce +1 0.026 697 599 60 23 
Ca 0 0.041 635 664 70 61 Cs 0 0.020 341 621 28 14 
ts aa 0.016 649 305 60 62 Ce fj 0.012 562 628 30 15 
Gs nd 0.005 000 062 00 63 1 =f 0.004 998 207 00 7 
C1 = 0.001 000 021 70 64 Co 8 0.001 000 000 00 1 
Co 4 0.000 100 000 00 65 1 =i 0.000 001 000 00 78 
Co 2 0.143 849 206 30 66 Cs +5 0.007 129 530 25 iT 


Pinboard Program: 


Starting point: u times 107 is in the accumulator. 
u should be shifted to meet this requirement. 


Oo W 67 
1, -R. 66 
2 = TELO 
ee. 

4 X 67 
Si 26k 
6 S: 15 
oe 03 
8:8 AGI3 
9". B 

10 W 68 
i” XX 68 
12 

13 

14 

15 


Result: 


e2" times 10-2 is in the accumulator; e" times 107? is 
in memory address 68. The results can then be shifted 
to meet the requirements of any problem. 


12 


Pinboard Program: 


Starting point: u times 107 is in the accumulator. 
u should be shifted to meet this requirement. 


0 W=s9 
1 R-79 
2 #H.10 
su 8 

4 X89 
5 + 7F 
6. S&S 7 
7 es 
8 Ws89 
9 B 

10 xX 89 
11 B 

124 RW 
13 == 189 
14 

15 


Result: 


e" times 10° is in memory address 89. The result 
can then be shifted to meet the requirements of any 
problem. 


15... 710" 
Range: 0<u<1 
Accuracy: Maximum error: +0.000 000 05 


Method: Reference—Hastings, page 144 
10% = [1 + cu + cu? + 2+ + + ca]? 


Number of Constants: 8 
Number of Temporary Memory Addresses: 1 
Switches Used: F (2 to 8) 


Time on E101: 7 seconds 


Constants: 
As Constant Assumed 

Constant Scale Factor Appears in Memory Memory Address 
Cr 0 0.000 932 642 67 90 
Ce 0 0.002 554 917 % 91 
C5 0 0.017 421 119 88 92 
C4 0 0.072 951 736 66 93 
C3 0 0.254 393 574 84 94 
Co 0 0.662 730 884 29 95 
C1 0 1.151 292 776 03 96 
1 0 1.000 000 000 00 97 


Pinboard Program: 


Starting point: u times 10° is in the accumulator. 
u should be shifted to meet this requirement. 


OA EE 2 
1 B 
Pe ¥)8) 
oe iol 
4 W 98 
5 "XK 198 
CES easels 
7 Sy 
8 U 04 
9 W 98 
10: 9 Av2t 
1 B 

PQ) xXeo8 
13 

14 

15 


Result: 


10" times 10-1 is in the accumulator. The result can 
then be shifted to meet the requirements of any problem. 


13 


16. log, u 


Range: 1 <u< 10 
Accuracy: Maximum error: +0.000 05 


Method: Reference—Hastings, page 126 


2 


— WV 10\ 2k 
log,)u = Lee Ss Conta ee 
k=0 


u+¥V10 


Number of Constants: 6 


Number of Temporary Memory Addresses: 2 


Switches Used: F (0 to 3) 


Time on E101: 5 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Cs | 0.025 432 750 00 50 
C3 | 0.027 738 390 00 60 
Ci —! 0.086 902 860 00 61 
1/2 —1 0.050 000 000 00 62 
V10 3 0.003 162 278 00 63 
rv 10 —3 0.006 324 556 00 64 


Pinboard Program: 


Starting point: u times 10-* is in the accumulator. 
u should be shifted to meet this requirement. 


On AGS 
1 B 

2.7) =x 164 
es)? 
4 R 52 
5 B 

6. 3052 
Lh OMSL 
8 H 10 
ad 3: 

10° XSF 
Bo =F IGE 
12 S 12 
13°. -W309 
14 

15 


Result: 


log, u times 10-1 is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 


17. 
Range: 1 << u< 10 


Accuracy: Maximum error: +0.000 05 


log, u 


Method: Reference—Hastings, page 126 


1 
log, u = M log,)u 


eg u — Y10\24 
“if $-C5] 
Number of Constants: 6 

Number of Temporary Memory Addresses: 2 
Switches Used: F (0 to 3) 


Time on E101: 5 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
Cs/M = 0.058 561 070 00 50 
C3/M 1 0.063 870 000 00 60 
C1/M 1 0.200 101 230 00 61 
1/2M =I 0.115 129 260 00 62 
V0 33 0.003 162 278 00 63 
2V 10 3 0.006 324 556 00 64 


Pinboard Program: 


The program is the same as that for log,, u. 


14 


18. 


Range: 1 <u< 10 
Accuracy: Maximum error: +0.000 000 003 
Method: Reference—Langdon, page 10 


log, u 


6 
log, u = log,c + ay Contr V7**1, where 
k=0 


moult 
pk usec 
Number of Constants: 10 


Number of Temporary Memory Addresses: 2 
Switches Used: F (0 to 6) 
Time on E101: 10 seconds 


Special Feature: increased accuracy over other log routine 
on this page. 


Constants: 
As Constant Assumed 

Constant Scale Factor Appears in Memory Memory Address 
Gis =2 0.000 572 283 27 80 
Ce 8 0.002 503 410 93 81 
& = 0.002 824 335 71 82 
Cs 2 0.004 001 930 33 83 
Cs = 0.006 666 617 10 84 
C1 =f 0.020 000 000 37 85 
loge ¢ of 0.011 512 925 46 86 
¢ 4 0.031 622 776 60 87 
Ic = 0.063 245 553 20 88 
Cis ay 0.004 105 970 44 89 


Pinboard Program: 


Starting point: u times 10-? is in the accumulator. 
u should be shifted to meet this requirement. 


Pinboard 1 Pinboard 2 
OF 87. 0 3B 

1 B 1 3X99 
2 — 88 2 + 86 
Bie Nee) 3 

4 R 99 4 

5. 1B 5 

6 «x99 6 

7 SB 7 

8 He 10 8 

9 R 89 9 

10 W 98 10 

11 x 98 ata 

12 + 8F 2 

13 Sets ED 

14 UW -0 10 14 

15 U 20 15 

Result: 


log, u times 10-? is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 


19. sinhu 


Range: —4.5 <u < 45 
Accuracy: Maximum error: +0.000 05 
Method: Reference—Taylor Series 
3 5 2n-1 
soi lela be saad ¢ 5 1b 
where n = 11 
Number of Constants: 11 
Number of Temporary Memory Addresses: 2 
Switches Used: F (0 to 9) 


Time on E101: 10 seconds 


20. coshu 


Range: —4.5 <u < 45 
Accuracy: Maximum error: +0.000 05 
Method: Reference—Taylor Series 


aohu=1t Ft tee toe, 
where n = 10 

Number of Constants: 11 

Number of Temporary Memory Addresses: 1 

Switches Used: F (1 to 10) 


Time on E101: 9 seconds 


Constants: Constants: 
As Constant Assumed As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address Constant Scale Factor Appears in Memory Memory Address 
1/21! 17 0.001 957 294 11 69 1 —4 0.000 100 000 00 69 
1/19! 15 0.008 220 635 25 70 1/20! 16 0.004 110 317 62 70 
1/11! 13 0.028 114 572 54 1 1/18! 14 0.015 619 206 97 11 
1/15! 11 0.076 471 637 32 72 1/16! 12 0.047 794 773 32 72 
1/13! 9 0.160 590 438 38 73 1/14! 10 0.114 707 455 98 73 
1/11! 7 0.250 521 083 85 14 1/12! 8 0.208 767 569 88 14 
1/9! 5 0.275 573 192 24 15 1/10! 6 0.275 573 192 24 15 
1/7! 3 0.198 412 698 41 76 1/8! A 0.248 015 873 02 16 
1/5! 1 0.083 333 333 33 11 1/6! 2 0.138 888 888 89 11 
1/3! —1 0.016 666 666 67 78 1/4! 0 0.041 666 666 67 78 
1 3 0.001 000 000 00 19 1/2! —2 0.005 000 000 00 19 


Pinboard Program: 


Starting point: u times 10-1 is in the accumulator. 
u should be shifted to meet this requirement. 


0 B 

1 W 68 
2 * 68 
5° eB: 

4 X 69 
5 Hedzo 
6 =Fin8 
7 W 67 
8 X=67 
9 ased8 
10: *W.06 
Lip Aaa 
12. -=5 

13. X68 
14 

15 


Result: 
sinh u times 10-* is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 


15 


Pinboard Program: 


Starting point: u times 1071 is in the accumulator. 


u should be shifted to meet this requirement. 


Oo B 

1 W 68 
2- “X68 
3» 4B 

4 X 70 
5- OT 
6— qe 7F 
7 W 68 
8 XxX 68 
2 S 19 
10 U 06 
11 + 69 
12 

13 

14 

15 


Result: 


coshu times 10-4 is in the accumulator. The result 
can then be shifted to meet the requirements of any 


problem. 


21. tanhu 


Range: —2 <u< 2 
Accuracy: Maximum error: 0.000 8 
Method: Reference—ElectroData 


c+ .245 


tanhu = where d = eee aca b, 


ee 

d= 1’ 
c= 'b-P .105; b =a? and a = "100 

Number of Constants: 3 

Number of Temporary Memory Addresses: 3 


Switches Used: none 


Time on E101: 1.9 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
105 0 0.105 000 000 00 00 
245 0 0.245 000 000 00 01 
100 0 0.100 000 000 00 02 


Pinboard Program: 


Starting point: u times 10 is in the accumulator. 
u should be shifted to meet this requirement. 


Oe B 

1 W 03 
2 eX, 03 
3 W 04 
4 + 00 
5.) 9-3 

eee wae cyl 
7 "05 
8 R 04 
Sine B 

LO X05 
aby + 02 
12 B 

13 R 03 
14 + 05 
15 


Result: 


tanh u times 10° is in memory address 05. The result 
can then be shifted to meet the requirements of any 
problem. 
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22. Vu 
(see note below) 
Range: Unrestricted, provided u times 10 is scaled so 


that n is an even integer and u times 10-* is less 
than 1. 


Accuracy: Maximum error in Y u times 10-*/? is less than 
0.000 000 000 01, if b = 6 in step 9. 


Method: Newton approximation: 


1 u 
Xin = A(x + =) 
i 


Number of Constants: 1 

Number of Temporary Memory Addresses: 2 

Switches Used: F (0 to 6) 

Time on E101: 1.17 seconds per iteration 

Special Feature: fewer program steps than routine on 


p. 17. 
Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
1/2 0 0.500 000 000 00 90 


Pinboard Program: 
Starting point: u times 10-" is in the accumulator. 
u should be shifted to meet this requirement. 


0” 38 

1 x 90 
2 W9i 
3 H 10 
ee il 
5 ==n9Z 
6 *-90 
J hee 
8a 

9: SS:lb 
10 U 04 
ala 

2 

13 

14 

15 

Result: 
Vu times 10-*/? is in the accumulator and in the B 
register. 


Note: u times 10-° must be scaled and placed in the ac- 
cumulator before control is transferred to the sub-routine. 
If the leading digit of u times 10 is the first or 
second digit after the machine decimal point, then 7 iter- 
ations (obtained by setting “b’ = 6 in step 9) will give 
the accuracy specified above. 

Study of the range of numbers whose square roots 
are to be taken will give some indication of the number 
of iterations required to achieve the desired accuracy; the 
“b” of step 9 can be pinned accordingly. 


23. Vu 


(see note below) 


Range: Unrestricted, provided u times 10™ is scaled so 
that n is an even integer and u times 10 is less 
than 1. 


Accuracy: Maximum error in Yu is 5 in the least signifi- 
cant digit of the accumulator. 


Method: Newton approximation: 


it u 
Xin = 2 (« =I a 
1 


Xiai = Vu when | x; — x;,,| < 0.000 000 000 05 
Number of Constants: 1 

Number of Temporary Memory Addresses: 3 

Switches Used: none 


Time on E101: 6 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
1/2 0 0.500 000 000 00 90 


Pinboard Program: 


Starting point: u times 10 is in the accumulator. 
u should be shifted to meet this requirement. 


Oo B 

je Te EC, 
2 W 91 
3 Ro 
4 + 92 
5 90 
G “er 92 
7 NOS 
Sir 4B 

Be RS 
10 —=— 592 
LU A 4 
12 A 24 
13 G03 
14> "R93 
15 


Result: 


Yu times 10-"/2 is in the accumulator, in the B 
register, and in memory address 93. 
Note: u times 10° must be scaled and placed in the ac- 
cumulator before control is transferred to the subroutine. 
If the leading digit of u times 10 is the first or 
second digit after the machine decimal point, then a good 
first approximation to ‘Vu times 10-*/? is obtained, and 
consequently the sequence of approximations converges 
very rapidly. 


17 


24. Multiplication of complex numbers 


Range: u + vi = (a + bi)(c + di) 
—.9999 <a,b,¢,d < .9999 


Accuracy: Depends on the accuracy of a, b, c, and d aad 
on the scale factors at which these numbers are 
entered. 

Method: (a + bi)(c + di) = (ac— bd) + (ad + bc)i 

Number of Constants: 4 

Number of Temporary Memory Addresses: 2 


Switches Used: none 


Time on E101: 2 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
a Both a and b must be entered with the 20 
b same scale factor; similarly, ¢ and d 21 
c must be entered with the same scale 22 
d factor. 23 


Pinboard Program: 


Starting point: a, b,c, andd in the assumed memory 
addresses. 


0 RY 21 
1 B 

2 R25 
3 W 24 
4 XX 22 
a. We25 
6 R 20 
7 B 

BO 5 
Oi Sie 25 
LO’ “W220 
11 X22 
12 — 24 
13: We2T 
14 

5 


Result: 


The real part of the result, (ac — bd), is stored in 
memory address 21 and in the accumulator. The imag- 
inary part of the result, (ad + bc), is stored in memory 
address 20. 


25. Division of complex numbers 


Range: u + vi = (a + bi)/(c + di) 
— 9999 <a, b, c,d < .9999 


Accuracy: Depends on the accuracy of a, b, c, and d and 
on the scale factors at which these numbers are 


entered. 
a+bi act+tbd , be—ad. 
Method: re oe ao pt 


Number of Constants: 4 
Number of Temporary Memory Addresses: 2 
Switches Used: none 


Time on E101: 4.5 seconds 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
a Both a and b must be entered with the 20 
b same scale factor; similarly, ¢ and d 21 
C must be entered with the same scale 22 
—d factor. 23 


Pinboard Program: 


Starting point: a, b, c, and —d in assumed memory 
addresses. 


Pinboard 1 Pinboard 2 
0 Re 21 Ka? 
1 B B 

2 25 x 22 
3 W 24 W 24 
4 De soz R 23 
Ba 25. B 

6 R 20 x23 
vs B + 24 
So ex 25 B 

9 + 25 R 20 
10 W 20 + 20 
LE x 22 Re 21 
52 — 24 eq ry)! 
13 We 21 

tA 20 

15 


Result: 


The real part of the result, (ac + bd)/(c? + d?), is 
stored in memory address 21. The imaginary part, 
(be — ad) /(c? + d?), is stored in memory address 20. 


C. Techniques Which Increase the Utility of the E101 
1. Example of Constant Increments Approach 
sin ku, and cos ku, for unit increments of k 
2. Example of Compound Functions Approach 
e" sinu 


3. Example of Multi-purpose Approach 


1. Example of Constant Increments Approach 
sin ku, and cos ku, for unit increments of k 
Range: All u, and k within the limits of the E101. 
Accuracy: Depends on the number of increments used. An 
example of the maximum error to be expected 
is as follows: uy = .04 radians; number of in- 
crements = 160; maximum error of +2 X 10-7, 
assuming sin Uy and cos uy are accurate to 10-1°. 
Method: 
sin kuy = cos uy sin (k — 1)u, 
+ sin uy cos (k— 1)uy 
cos kuy = cos uy cos (k — 1)u, 
— sin Uy sin (k — 1)uy 
Number of Constants: 4 
Number of Temporary Memory Addresses: 2 
Switches Used: none 
Time on E101: 2 seconds for each k 


Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
sino 0 0.000 000 000 00 40 
sin Uo 0 sin Uo times 10° 4] 
COS Uo 0 COS Uo times 10° 42 
cos 0 0 1.000 000 000 00 43 


Pinboard Program: 
Starting point: sin 0, sinuy, cos0, and cos u, times 
10° are in assumed memory addresses. 


10) R 41 
iD 

2 X 43 
3 W 44 
4 X 40 
5 W 45 
6 R 42 
7 <3 

8 xX 40 
9 +44 
10 W 40 
11 X 43 
12 — 45 
13 W 43 
14 

lb, 


Result: 
sin ku, times 10° appears in memory address 40; 
cos ku, times 10°, in memory address 43. 


2. Example of Compound Functions Approach 
e" sinu 
Range: —2 <u< 2 
Accuracy: Maximum error: +0.000 1 
Method: Reference—Taylor Series 


e'sinu =u + u? + =e Out — <u 
a au = ut + Ou8 + iru 
+ = ye a a yit 
Number of Constants: 11 
Number of Temporary Memory Addresses: 1 
Switches Used: F (0 to 10) 
Time on E101: 10 seconds 
Constants: 
As Constant Assumed 
Constant Scale Factor Appears in Memory Memory Address 
32/10! 5 0.881 834 215 04 30 
16/9! 4 0.440 917 107 52 31 
0 0 0.000 000 000 00 32 
—8/7! 2 0.158 730 158 72— 33 
—8/6! 1 0.111 111 111 «11— 34 
—4/5! 0 0.033 333 333 32— 35 
0 0 0.000 000 000 00 36 
2/3! —2 0.003 333 333 33 37 
1 3 0.001 000 000 00 38 
1 —4 0.000 100 000 00 39 
32/11! 6 0.801 667 468 48 49 


Pinboard Program: 
Starting point: u times 10+ is in the accumulator. 
u should be shifted to meet this requirement. 


Oo B 

T 3x49 
2 Et 10 
bie ape ele! 
4 W 48 
5 xX 48 
6. (Sito 
eo AUEOS 
8 

9 

10 

11 

12 

13 

14 

15 


Result: 
e" sin u times 10-5 is in the accumulator. The result 
can then be shifted to meet the requirements of any 
problem. 
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3. Example of Multi-Purpose Approach 


This program, initially mentioned on page 2 of Section I, 
is an example of a multi-purpose subroutine that can be 
used to compute sin u, cos u, arctan u, and e~”. It is assumed 
here that these four functions occur in one problem. The 
illustration shows how to arrange the memory and how to 
transfer into, and out of, the multi-purpose subroutine. 

Further, the illustration assumes that sin u is to be called 
for in pinboard 1; cos u in pinboard 2; arctan u in pinboard 
4; and e™ in pinboard 6. 


Constants: 


The constants for sin u are stored in row 1 of the mem- 
ory; the constants for cosu are stored in row 2 of the 
memory; the constants for arctan u are stored in row 4 of 
the memory; and the constants for e are stored in row 6 
of the memory. It is important to note that the number of 
the row in which the constants for a particular function are 
stored is the same as the number of the pinboard in which 
the function is called for. 


MEMORY* 
Di ee ee eS ae ee 
0 
HMA & Aft 1A! + SAietetihes 
ee Wit = b 1H, 8 Hi 6. Y 
3 
ee ee a Re OE ES oe See 
5 
6 


See es OE Os PM dees Segue ea 


*For range, accuracy, scale factor, time on E101, and con- 
stants as they appear in memory, of sin u, see page 6; of 
cos u, see pages 7 and 8; of arctan u, see pages 10 and 11; 
and of e", see page 12. 


Program: 


The program for the multi-purpose subroutine is shown 
in pinboard 8. We assume, of course, that each time we 
enter the subroutine, u has been properly scaled and placed 
in the accumulator. 

Just as in any multiple use of a subroutine, the E switch 
is homed to the number of the pinboard to which the 
control is to return when the subroutine has ended. But, 
because of the arrangement of the memory, as designated 
above, the E switch is also used to reference the proper 
row of constants. The F switch is then used to iterate along 
the chosen row. 


Pinboard No. 


U80 


9 (sin u) 


This program is an illustration of a technique. Other polynomial approximations may, 
of course, be adapted to the same general procedure. However, one word of advice may 
be in order: it is perhaps most convenient first to arrange the memory and then adapt the 
program to exit from the proper pinboards, rather than to arrange the memory to corte- 


spond to the program exits. 
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